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Abstract 

Massive gravity previously constructed as the spin-2 quantum gauge 
tlieory leads in the mass zero limit to a modification of general relativ- 
ity. As a relic from the massive theory a vector field survives which 
couples to the metric only. However, the coupling does not involve 
covariant derivatives so that must be considered as 4 scalar fields. 
We analyze the static, spherically symmetric solutions of this theory. 
From the corresponding geodesies we find the circular velocity profile. 
Interpreting this as coming from a dark density profile, the theory pre- 
dicts a flat density core for r ^ 0. But the dark density profile is not 
universal. 
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1 Introduction 



For some years we have studied massive gravity on the basis of a cohomolog- 
ical formulation of gauge invariance in terms of asymptotic free fields. This 
method works for massless and massive gauge theories equally well [1], the 
Higgs mechanism is not used. The massive spin-2 theory was first investi- 
gated in [2]. The most elegant way to obtain the theory is by assuming the 
gauge invariance condition for all chronological products in the form of the 
descent equations [3] . These give the total interaction Lagrangian including 
ghost couplings and the necessary coupling to a vector-graviton field v^. 
This vector-graviton field is of fundamental importance because without it 
a massive spin-2 gauge theory is impossible. So the main purpose of this 
paper is to analyze the physical consequences of this new actor in the drama 
of gravity 

Einstein's theory describes reality very well on the scale of the solar 
system, that means on distances up to 100 AU. On the scale of galaxies, i.e. 
kilo-par sec, one observes deviations from general relativity which are usually 
ascribed to hypothetical dark matter. In the tensor-four-scalars theory they 
come from the vector-graviton field v'^ which must be considered as four 
scalar fields in the classical theory. Since the scalar fields v'^ and gravity 
are coupled by derivative couplings it is the spatial variation of which 
matters. Prom the two length scales above it is plausible that the mean 
spatial variation of v'^ over a galaxy is 10^ times bigger than in the solar 
system. This explains why the deviations from general relativity are only 
observable on the galactic scale or bigger. 

The paper is organized as follows. In the next section we recall the origin 
of the modified gravity theory. In section 3 we introduce a galaxy model and 
construct the static, spherically symmetric (vacuum) solution. The latter 
was already tried in a previous paper [4] with an ansatz containing two free 
metric functions. This was not good enough, we need three functions. In 
section 4 we compute the power series expansions of the solution for r — >■ 
and r — 7- oo. In section 5 we consider the geodesies corresponding to our 
metric and specialize to circular motion. This gives an expression of the 
circular velocity in terms of the metric functions. Interpreting this as the 
velocity curve of dark "matter" we find a flat density profile for r — >■ 0. In 
section 6 wc discuss the degeneracy in the field equations. As a consequence 
we find that the dark density profile is not universal. 
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2 Massive quantum gravity and the tensor-four- 
scalars theory 

Our starting point is the massive spin-2 quantum gauge theory on Minkowski 
space which we have called massive gravity for short. The basic free asymp- 
totic fields are the symmetric tensor field h^^^{x) with arbitrary trace, the 
fermionic ghost u^^{x) and anti-ghost u^^{x) fields and the vector-graviton 
field v^{x). They all satisfy the Klein-Gordon equation 

(□ + m'^)hi^'' = = {D + m'^)ui^ = (□ + m'^)u'' = {D + rr?)v^ (2.1) 

and are quantized as follows [2] [3] 

[h'^^{x),h^''{y)] = -l(ry"/^r/^^ + ry^'^r?'^^ - <^r?'^'^)L>^(x - y) (2.2) 

K(x),«^(y)} = i77^^Z^„(x-2/) 
K(x), v^{y)]='-rii'^Dm{x-y), 

and zero otherwise. Here, D^{x) is the Jordan-Pauli distribution with mass 
m and r]'^^ = diag(l, —1, —1, —1) the Minkowski tensor 

The gauge structure on these fields is defined through a nilpotent gauge 
charge operator Q satisfying 

Q2 = 0, Qn = {) (2.3) 
where Q, is the Fock vacuum and 

dg/jM^ = [g, /i/^^] = -^-{d^u^ + di'u'' - r?^''5«u") (2.4) 

dQU^ = {Q,u} = (} 

dQU^' = {Q, u^'} = {{d^h^" - mv") (2.5) 

dQV^'=^ [Q,v^'] = -^mu^". (2.6) 

The vector-graviton field is necessary for nilpotency of Q. In fact, in 
order to get dqU = from (2.5), the additional term —mv^^ is indispensable. 

The coupling T{x) between these fields follows from the gauge invariance 
condition [2] 

dQT{x) = id^T'^ix) (2.7) 
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where T and are normally ordered polynomials with ghost number and 
1, respectively. In addition we may require the descent equations 

SqT" = [Q, T"] = idfiT''^ (2.8) 

[Q,T'^^]=id^T'^^^ (2.9) 

where the new T's are antisymmetric in the Lorentz indices. The essentially 
unique coupling derived from (2.7-9) is given by [3] 

-h'^^h^ph+^h^). (2.10) 

Here h = h'^ is the trace and a coupling constant is arbitrary. The quartic 
couplings follow from second order gauge invar iance and so on. 

We consider the limit m — t- in the following. The massless limit of 
massive gravity is certainly a possible alternative to general relativity. The 
new physics comes from the surviving coupling term of the vector-graviton 

n = AK'^daVxdpv^. (2.11) 

To be able to do non-perturbative calculations wc look for the classical 
theory corresponding to the coupling (2.10). It was shown in [1] sect. 5. 5 
that the pure graviton couplings hdhdh correspond to the Einstein-Hilbert 
Lagrangian 

2 , 

-^EH = — nV-gR 

in the following sense. We write the metric tensor as 

and expand Leh in powers of k. Then the quadratic terms 0{k^) give the 
free theory, the cubic terms 0{k^) agree with the pure graviton coupling 
terms in (2.10) up to a factor and and quartic and higher couplings follow 
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from higher order gauge invariance, so that the full non-linear structure of 
general relativity is recovered. 

To obtain (2.11) we must add 



to Leh- One may be tempted to write covariant derivatives Vq instead 
of partial derivatives in order to get a true scalar under general coordi- 
nate transformations. But this would produce quartic couplings containing 
v\ and such terms arc absent in the quantum theory ([2], eq.(4.12)). For 
the same reason the Lorcntz index A in vx is raised and lowered with the 
Minkowski tensor rjn^, but all other indices with ga^- Both together means 
that the vector graviton field v\ should he considered as four scalar fields 
in the classical theory. Then (2.12) is a scalar under general coordinate 
transformations as it must be. 

Of course we want to include ordinary matter in the theory which we 
describe by a complex scalar field (p of mass M. In the Klein-Gordon equa- 
tion for this field we retain Planck's constant U different from 1. In the 
end we want consider the classical limit ti. 0. Then in the semi-classical 
approximation the solution of the Klein- Gordon equation is of the form 



and it has the following interpretation: q(x) is the density of particles in 
space and time and 

is the momentum of the particle at x. S{x) is the Hamilton-Jacobi principal 
function. Nobody forbids to take for M the solar mass so that the particles 
are the stars in a galaxy, In this way we get a simple model of a galaxy. It 
is indeed simple because all stars have the same mass as the sun and gas 
etc. is neglected. On the other hand the coupling of the complex scalar 
field to gravity is completely known from the analysis of gauge invariance 
[1]; we do not need equations of state as in more realistic models. The total 
Lagrangian of our model now reads 




(2.12) 




(2.13) 




(2.15) 
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The Euler-Lagrange equations for the Lagrangian (2.15) give the system 
of coupled field equations. Variation of g^'^ yields the modified Einstein 
equations 



^2 



-d^vxd.v^ + ^9,^,9''^ davxdpv^} (2.16) 

The variational derivative with respect to v^^ gives the wave equation in the 
metric 9°^^ 

^daiV^g^'^d^v^) = 0. (2.17) 

Finally, the variation of where is not varied gives the Klein-Gordon 
equation in the metric 9°^^ 



da{V^9"'^dfsip) + M\ = 0. (2.18) 



From the point of view of quantum gauge invariance this rather simple 
tensor-four-scalars theory has the same right for being considered as funda- 
mental as Einstein's theory. The latter is the somewhat exceptional m = 
theory; the former is the massless limit of the massive spin-2 gauge theory. 



3 Static spherically symmetric solutions 

As a first step to analyze the tensor-four-scalars theory we study vacuum 
solutions of the field equations (2.16) (2.17) neglecting the normal matter 
((/? = 0). Of course the Schwarzschild solution of ordinary general relativity 
is also a solution of our theory with = const. This solution describes 
reality very well in the solar system. But on larger scales the spatial variation 
of becomes observable, so that the source terms depending on derivatives 
of must be taken into account. In this sense, Einstein's theory can he 
considered as the local version of the tensor-four-scalars theory. It is our 
aim now to study solutions with non-constant v^. We choose spherical 
coordinates 

= ct, = r, x^ = x^ = (p (3-1) 
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and assume the metric ds'^ = g^ydx'^dx'^ to be of the fohowing static spher- 
ically symmetric form 

ds^ = e'^c^dt^ - e^dr'^ - r^e^id-d^ + sin^ Mcp^), (3.2) 

where a(r) , b{r) , c(r) are functions of r only. Then the determinant of the 
metric is equal to 

g = -e^"^ sin^ ^, d = a + b + 2c (3.3) 
and we obtain the following non-vanishing Christoffel symbols: 

r?o = ^, rlo = p-' (3.4) 

r}i = |, r^2 = -(r + ^c')e-'' (3.5) 

ri3 = sin2^r^2, r?2 = i + ^ = r?3 (3.6) 

rig = - sin cos ^, rig = cot (3.7) 

The prime always means partial derivative with respect to r. 

The next step is the computation of the Ricci tensor. Only the diagonal 
elements are different from zero: 

i?oo = ^e^-\a" + - ^a'b' + a'c' + ^a') (3.8) 

i?n = -hia" + 2c") + ^(a' + 2c' + J) - ^ - y - ^c' (3.9) 

i?22 = e'^-^i-l - ^c" - r(2c' + ^^) - jc'{a' - b' + 2c')] + 1 (3.10) 

i?33 = sin2?9i?22. (3.11) 
This gives the following scalar curvature 

R = g'^^'R^,,, 

= e-'' [a" + 2c" + ^ - ^ + a'c' + -a' - -6' - b'c' + ^c'^ + -c'] + - (e'^ - e"-) . 
2 2 r r 2 r r 

(3.12) 
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Now we arc ready to calculate the left-hand side of the modified Einstein 
equation (2.16) 

Gfj,i, = Rixv — Ig^i/R. (3.13) 
It is convenient to raise one index by multiplication with g^'^: 



G'o = e-\-c" - + ^ + ^ + (3.14) 

4 z J. j.^ j.^ 



+ ^ (3.15) 

«i ^ -'(444.^-^.^4.^^) ^ oi. ,3.1.) 

Before we calculate the right-hand side of the modified Einstein's equa- 
tion we consider the wave operator in the metric (2.17) neglecting the time 
derivatives Oq: 

pd/2-c 

'''' (3.17) 

Then the field equation for a static becomes 

9/d/2-b^2d\^ = ed/^-OLW (3.18) 

where is the quantum mechanical angular momentum operator squared. 
Since we restrict ourselves to spherically symmetric solutions the right-hand 
side vanishes. Then only the time component = {vo{r), 0,0,0) can be 
different from 0, otherwise we get non-diagonal elements in the metric. The 
equation (3.18) can now be integrated once 

^o(0 = ^e^-'/', (3.19) 

where A is an integration constant. With this result the modified Einstein 
equations (2.16) assume the following form: 



2 r4 



GO = -^e-'^fl^e'"'-'' (3.20) 



= |^e^-<^ (3.21) 
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.2 _ ,9_:^^ „b-d 

2 r 

where 



G'l = -^TXe'-' (3.22) 



9 = (3-23) 

is essentially Newton's constant. Using the above results for the Einstein 
tensor and multiplying by we get the following three differential equa- 
tions 

c" = + h'c' + -{b'- 3c') + l(e^-^ - 1) + ^^e""'-^ (3.24) 

4 2 r 2r^ 



= '-a'c' + l(a' + cO + ^ + ^(l - e^"^)+^,e^^-' (3-25) 

II , It 1 /W / o '\ ■'■ / '2 /i/ I / / I,' / 1 I2\ I 26— d /o r,c\ 

a +c = -{b — a — 2c ) — - {a — ab+ac— bc+c ) + —^e , (3.26) 
where 

a = gA^ (3.27) 
is another form of the constant of integration for VQ{r). 



4 Power series expansions 

At first sight the three equations (4.24-26) for the three functions a{r),b{r) 
and c(r) seem to look quite awful because of the r~^-singularity which comes 
from the t;o-function. To study the solution for small r we try to solve the 
equations by power series 

CXD 

air) = J2 (^nr'' 4.1) 

n=0 

OD 

b{r) = E (4-2) 

n=0 

However in c(r) we need a logarithmic term 

oo 

c(r) = -2 log- + Ecnr-" (4.3) 

'^^ n=0 
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in order to cancel the singularities. We expand everything in powers of r, 
for example 



(7) ^e*""""" [1 + (^1 - ci)r + (62 - C2 + 



bl + cl 



61 Cl) 



+ . 



(4.4) 



r \ 4 



-a-2c ^ 



-ao—2cQ 



1 + (61 - ai - 2ci)r+ 



al + hl + 44 



— aibi — 2b\Ci + 2aiCij 



For zeroth order we introduce the two parameters 

1 



060— CO 



It is convenient to subtract (3.24) from (3.26) 

a" = ^ - ^(a' - 6' + c') + -{d -a') + \{l- e^"^) 
4 2^ r ' 

^2r4'' 



(4.5) 

(4.6) 
(4.7) 



(4.8) 



and to use this equation instead of (3.26). Then to lowest order 0{r~^) all 
three equations (6.2.24), (6.2.25) and (4.8) are identically satisfied without 
giving a restriction on ai,6i,ci. In next order 0(r°) we obtain 



02 



ai = -| + -(2/3-7) 

c| _ of aihi _ aici _ P_ 1_ 
8 4 4 4 2 4 



(4.9) 



16 



61 Cl 



+ (2/3-7) 



61 (2/? -7)' P , 7 



3 2 . ^ici , . 7 

C2 = -qCi H h — + T- 



4ci 4cf 
2 4 



(4.10) 
(4.11) 



In the order 0{r) we find 



= -32^^^ - - 1^ - 24^^^ - 16^^ + 12^ - 24^ + 6^ - 2^+ 
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5 7^_^6i^5i_7^6i 2^_^ /^T^ 7' 



24 ci 2 cf 2 cf 8 c? 3 cf cf 2cf 12cf 

.^61 + ^61-^ + ^^-^^. (4.12) 
4^8 12 3 ci 6 ci ^ ^ 



and 



C3 = l6^^^-^^^^^+^+^^^^-12^'^^-4^^-6^+12^+4^^+8^^- ^^''^^ 

However, the equation (3.25) without second derivatives is identically sat- 
isfied and gives no further relation. The same is true in O(r^), O(r^) and 

O(r^) (see [1]). 

Assuming that this property holds in all orders, then all n > 1, c„, n > 
2 are determined by ci and 6„, n > 1. That means the function b{r) remains 
completely free. We get a class of solutions with one arbitrary function b{r). 
We shall return to this point in section 6. 

For r — )• 00 we expect the solution to be asymptotically flat because it 
should represent an isolated model galaxy. Therefore we set up an expansion 
of the form 

« = E ^ (4-14) 

n=l 

^=E^ (4-15) 

n=l 

c=j:% (4.16) 

n=l 

Then the exponentials are expanded as follows 



e- = l + ^ + l(5.-C. + l(5.-C.)^) + 



+l(i?3 - C3 + (Si - Ci)(S2 - C2) + ^(5i - Ci)3)+0(r-4) 
= 1 + " " + ^^(^2 - ^2 - 2C2 + \{B, -A^- 2C,f) + 
+^{B3 -A3- 2C3 + (Bi -Ai- 2Ci){B2 -A2- 2C2)+ 
+hB^-Ai-2Cif)+0{r-^) 
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Substituting this into (3.24), (3.25) and (4.8) we get from 0(r"^) 

Ai = -Bi. (4.17) 
In order 0(r~^) it follows from (3.24) 



C, = :?i-Ml_B,-^2 + « (4.18) 

2 2 4 2 ^ ^ 



and (4.8) determines 



A, = -^ + ^. (4.19) 



Prom (3.25) we obtain the quadratic equation 

= (Si - Ci)2 + C1-B1 (4.20) 

hence 

Ci = Bi, or Ci = Si - 1 (4.21) 

so that we find two possible asymptotic solutions. 
In next order 0{r~^) (3.24) gives 

BjCi B,B2 B,Cl B, 

B^ B, Cl_a (4.22) 
2 2 8 2 24 ^ ^ 

and from (4.8) we get 

A Bf B1B2 BiCf 

^ 18 6 54 2 6 

+ + + (4.23) 

27 6 108 36 54 108 ^ ^ 

Then the equation (3.25) is identically satisfied for both solutions (4.21). 
The same is true in 0{r^^) and 0(r^^) (see [1]): from (3.24) one finds C„, 
An is given by (4.8), and (3.25) gives no further restriction on the solution. 
As in the expansion for small r, b{r) is completely free and determines a(r) 
and c(r); but there is a small difference: ci remains free for small r, whereas 
Ci is fixed according to (4.21). So if one tries to match the two expansions 
by numerical integration of the field equations one better starts from large 
r and integrates backwards to small r. Counting the degrees of freedom it 
should then be possible to match the solution. To learn how this can be done 
in a realistic case we discuss the connection of the metric with observable 
quantities in the next section. 
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5 Geodesies and rotation curves 



As in ordinary general relativity wc assume that test particles like gas or 
stars in a galaxy move along geodesies. The geodesic equation reads 

where s is an affine parameter. We restrict ourselves to motion in the plane 
•d = 7r/2. Using the Christoffel symbols (3.4-7) we have the following three 
equations 

cFct idctdr ^ ('5 2) 

ds"^ ds ds 

We further restrict ourselves to circular motion r =const, so that dr/ds = 0. 
Then from (5.3) we get 



2 \ ds J V 2 / ydsJ 

This implies 

(5.6) 



/ dip \ a' 



\dct' 2r + r^c' 

The left-hand side in (5.6) is essentially the rotation velocity v or angular 
velocity u 

^ = ^ = ^, (5.7) 

act c cr 

so that we obtain the important result 

-e"-^. (5.8) 



2 / 

V ar 



c 



:2 2 + c'r 



We observe with satisfaction that the measurable quantity v depends on the 
metric functions a(r) and c(r), only, not on b{r) which was arbitrary in the 
power series expansions. 
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Since we have neglected ordinary matter in our model we must interpret 
v{r) (5.8) as the rotational velocity of the dark matter halo. Usually this 
halo is described by a density profile g{r) or mass profile 

r 

M{r) = 47r ^(r)r^ Sx. (5.9) 



Assuming Kepler's third law 

V- = ^ (5.10) 

r 

the dark matter density profile is given by 

This profile is usually determined from numerical simulations on the basis of 
Newtonian dynamics ([5] and references given therein). On the other hand 
the total density (dark plus ordinary matter) can be measured by means of 
the rotation curve v{r) of the galaxy. We now have a third approach to this 
quantity with our tensor-four-scalars theory. 

Of particular interest is the behavior of Q{r) for small r because it is 
difficult to determine this by numerical simulation. Substituting the power 
series expansions (4.1) and (4.3) 

ra' = air + 2a2r^ + Saar^ + . . . 

rd = -2 + cir + 2c2r^ + Scar^ + . . . (5.12) 
into (5.8), we get the following expansion for the circular velocity profile 



v"^ _ 2 + 2a2r + 3a3r2 + . . . e"0"'=o 
(? ci + 2c2r + 3c3r2 + . . . r 



.2 _|_ 0^„». I 'i.n„ri -X- ^ ^ ' ^ ^ ' 



Since the 2 in the denominator in (5.8) is cancelled by the logarithmic terms 
in (4.3), the circular velocity v has a linear behavior v{r) r for r — > 0. 
This can be traced back to the action of the vector graviton or scalar field in 
(3.19). For the density profile (5.11) this means that the tensor-four-scalars 
theory predicts a flat dark density core Q{r) ~ r^ for r — t- 0. On the other 
hand the older numerical simulations of Navarro, Frenk and White [6] have 
given an inner cusp as steep as r"^. 
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In this situation one hopes that observations can distinguish between the 
different predictions. Probably the best galaxy for this purpose is M33. This 
is a low-luminosity spiral galaxy in the Local Group which is dark-mattcr- 
dominated. A detailed study of the radial distribution of visible and dark 
matter was made by E. Corbelli [7]. Corbelli has fitted the measured data 
using dark density profiles with both fiat and cuspy cores. Unfortunately, 
the fits were equally good. But in later work by Donato et al.[8] the fiat 
dark density core was confirmed. 



6 Reduction of the equations and discussion 

In section 4 we have found that the metric fvmction 6(r) was not determined 
by the field equations (3.24) (3.25) (4.8). Another way to see this is the 
following. We solve equ.(3.25) for b{r) 

6 = c + 2 log r - log(l - ^e"""') + 



a c 



a' + c' c'2 1 



+ + ^ + - + (6-1) 

This enables us to eliminate h in the remaining two equations (3.24) and 
(4.8). Then from both equations we arrive at the same single equation 
between a{r) and h{r) 

a' \ r' zr^ — a exp(— a — cj 

This shows the degeneracy in the system of vacuum field equations. 

A similar degeneracy is known from the Schwarzschild solution with the 
difference that there the two remaining equations allow to calculate the two 
metric functions. Including the dark sector by means of the -y-field, the 
vacuum problem without normal matter becomes undetermined. To maJce 
it determined we need a further equation for a and c. For this purpose we 
can use equ.(5.8) for the rotation velocity in the form 

a'e" = -^(2 + rc')e^ (6.3) 

Hence, if the rotation velocity is known the metric can be calculated. This 
implies that there is no universal dark density profile in contradiction to ear- 
lier numerical simulations [6]. One expects that the degeneracy in the field 
equations is lifted if the ordinary matter is included in the form described 
in section 2. But this is a much more difficult problem. 
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